Ab initio studies of the spin-transfer torque in tunnel junctions 
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We calculate the spin-transfer torque in Fe/MgO/Fe tunnel junctions and compare the results to 
those for all-metallic junctions. We show that the spin-transfer torque is interfacial in the ferro- 
magnetic layer to a greater degree than in all-metallic junctions. This result originates in the half 
metallic behavior of Fe for the Ai states at the Brillouin zone center; in contrast to all-metallic 
structures, dephasing does not play an important role. We further show that it is possible to get 
a component of the torque that is out of the plane of the magnetizations and that is linear in the 
bias. However, observation of such a torque requires highly ideal samples. In samples with typical 
interfacial roughness, the torque is similar to that in all-metallic multilayers, although for different 
reasons. 

PACS numbers: 73.63.-b,75.70.Cn,75.30.Et,71.15.Mb 
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The discovery of giant magnetoresistance (GMR) in 
spin- valve systems [J Q and the rediscovery of the tun- 
neling magnetoresistance (TMR) in tunnel junctions [Sj, 
Q has led to applications such as hard-disk read heads, 
sensors, and storage elements in magnetic random access 
memory (MRAM). Both GMR and TMR occur in junc- 
tions in which two ferromagnetic leads are separated by 
a spacer layer, which is a non-magnetic metal in case of a 
GMR junction and a tunnel barrier for a TMR junction. 
The resistance R and therefore the conductance g of these 
junctions are a function of the relative angle 6 between 
the magnetizations of the ferromagnetic leads. The mag- 
netoresistance (MR) ratio is [i?(180°) - R{0°)]/R{0°). 
Typical GMR ratios are in the range of 50 % 5] and 
can be explained by spin dependent interface scattering 
in a semiclassical picture [fl lZ, ,8i] . TMR ratios can exceed 
several hundred percent in crystalline Fe/MgO/Fe tunnel 
junctions [1, as predicted theoretically [ill. 

These high TMR values in crystaUine Fe/MgO/Fe 
junctions originate in the symmetry-dependent transmis- 
sion probabilities through the MgO barrier at the Bril- 
louin zone center (F point) combined with the exchange 
splitting of the Fe band structure d, [IH . States which 
have the full rotational symmetry of the interface, said 
to have Ai symmetry, decay the most slowly in MgO 
and hence dominate the tunneling current. At the Fermi 
level, Fe is a half metal at the F point for states with Ai 
symmetry, having only majority states. The half metallic 
nature of the states that dominate the tunneling leads to 
a much higher current for parallel than for antiparallel 
alignment of the magnetizations. 

Spin-transfer torque, an effect predicted by Slon- 
czewski 1^ 14 1 and Berger 15 1, can be used to switch 
the magnetic orientation of ferromagnetic layers in GMR 
and TMR devices. For this purpose a current is driven 
through the sample and becomes spin polarized in one 
ferromagnetic layer. This polarization persists going 
through the spacer layer and entering the other ferro- 



magnetic layer. If the spins of the polarized current are 
not aligned with the magnetization, they precess around 
it. This precession in turn creates a torque on the magne- 
tization and can reverse the magnetization if the current 
is high enough. There is currently significant interest in 
understanding spin-transfer torques in tunnel junctions 
as a way to allow the development of MRAM applica- 
tions |16j . 

Spin transfer torques are well understood in all- 
metallic trilayer structures [17]. There, the current is 
carried by electrons over the whole Fermi surface. In fer- 
romagnetic layers and at their interfaces, electrons pre- 
cess at different rates and the components of the spins 
transverse to the magnetization rapidly become out of 
phase from each other. Two properties of the torque fol- 
low from the strong dephasing of the electron spins [18i] . 
First, the spin transfer torque largely occurs at the in- 
terfaces between ferromagnetic layers and non-magnetic 
layers. Second, it is largely in the plane defined by the 
magnetizations of the two ferromagnetic layers. These 
properties are used in almost all modelling of dynamics 
of GMR devices [17] and have been used without change 
in the modelling of TMR devices. 

In typical tunnel junctions, the current and spin cur- 
rent are carried by a small fraction of the Fermi surface 
and dephasing is greatly reduced. Here, we use ab ini- 
tio calculations to compute the spin transfer torques in 
Fe/MgO/Fe tunnel junctions. We show that in spite of 
the reduced dephasing, the torque is still approximately 
confined to the interface and the in-plane torque is much 
larger than the out-of-plane torque (see Fig.[T]). 

We treat the same structure as in previous studies on 
Fe/MgO/Fe [3] taking into account the experimentally 
observed relaxation of the Fe layer next to the inter- 
face [2^. The junctions consist of an Fe fixed layer, an 
MgO barrier, and an Fe free layer embedded between 
semi-infinite Cu in a bcc-Fe structure. The potentials 
are calculated self-consistently using a screened Korringa 
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FIG. 1: (Color online) Top: Schematic geometry. The left 
Ml and right Mr magnetizations of the tunnel junction are 
taken to lie in in the xz plane, at a relative angle 9 (here 
taken to be 90 °). The spin transfer torque acts perpendicular 
to each magnetization and can be divided into the in-plane 
torque ry, which lies in the xz plane, and the out-of-plane 
torque t± , which points into the y direction perpendicular to 
the plane defined by Ml and Mr. 



Kohn Rostoker multiple scattering Green's function ap- 
proach. We obtain potentials for junctions with differ- 
ent free layer thicknesses using the frozen potential ap- 
proximation and the potentials from a calculation with 
20 Fe monolayers in each magnetic layer. We calculate 
the linear-response torque using the magnetic moment 
S'irL{Ep) of those electrons that contribute to the cur- 
rent. The torkance dr/dV is the variation of the torque 
r with the voltage V . The torkance acting on layer i 
is 
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where A is the exchange field pointing in the direction 
of the magnetization of the corresponding layer, A is the 
area of the in-plane unit cell, and go = /h is the quan- 
tum of conductance. For a description of our implemen- 
tation of the no n-eq uilibrium Green's function technique 
see Refs. [HI l23l The total non-equilibrium magneti- 
zation 5mi{Ep) = (l/2)[(5m^(i?F) - 6rnf{EF)] contains 
separate contributions from electrons incident from the 
left and holes incident from right (for a positive bias). 

Fig. [2] shows our ab initio calculations of the torkance 
as a function of the relative angle 9 between Ml and Mji 
for different thicknesses of the free layer. Both the com- 
ponent in the plane of the two magnetizations, dr^^/dV, 
and the component out of the plane dr^/dV are al- 
most perfectly sinusoidal. For current biased applica- 
tions, the torque per current df/dl is of greater inter- 
est. This is related to the torkance by the conductivity 
g — dl /dV, g{df/dl) = dr/dV. In tunnel junctions with 
very high TMR ratios the conductivity depends strongly 
on the angle between the magnetizations so that df / dl is 
highly asymmetric. The critical voltages and currents for 
switching out of the parallel and antiparallel states are 
proportional to the inverses of the slopes of these curves 
at0 = O° and 180° respectively. 

The bottom left panel of Fig. [2] shows that both the 
in-plane and out-of-plane torkance oscillate as a func- 
tion of the Fe free layer thickness. The oscillations fit a 
sine curve with a period that is incommensurate with the 
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FIG. 2: (Color online) Top: Torkance (left) and torque per 
current (right) acting on the free layer as a function of the 
relative angle 6 between the magnetizations of the ferromag- 
netic electrodes. Results are shown for different thicknesses of 
the Fe free layer as indicated in the figure. Bottom: Torkance 
(left) and torque per current (right) as a function of the free 
layer thickness for an angle of 90 °. All calculations have a 
20 monolayer fixed layer and a 6 monolayer barrier. 



lattice spacing (hence the apparent beating of the ampli- 
tude). This period is very close to the Fermi wave vector 
of the Ai band in Fe at the F point. This agreement 
of the periods indicates that the important states are 
located close to the Brillouin zone center. The conduc- 
tance shows similar oscillations so that the torque 
per current in the right bottom panel of Fig. [2] is largely 
independent of the thickness. However, there is a phase 
shift between the oscillations in the in-plane and out-of- 
plane torkance so that the oscillations in dT± /dl are even 
stronger than they are in the torkance. 

Similar oscillations of the conductance and magnetore- 
sistance have been observed as a function of the thick- 
ness of an additional nonmagnetic layer inserted next to 
the barrier [Hm. For Fe/MgO/Fe, the origin of these 
oscillations is subtle. As we show below, the torque is 
largely restricted to the interface because the minority 
component is evanescent. However, the torque changes 
as a function of the free layer thickness due to coherent 
effects: specifically, phase shifts in the majority channel 
due to quantum well states. 

Note that in Fig. [5] the torque per current is very close 
to fiB/e- when the magnetizations are at an angle of 90 °. 
This is the expected behavior for a simple model of a 
junction between two half metals, for which the polar- 
ization of the current is 100 %. Each electron spin that 
traverses the junction rotates by 90°. This change in an- 
gular momentum is shared by the fixed and free layers. 
Since the torque must be perpendicular to the magne- 
tizations, the change in angular momentum supplied to 
each layer must be hb- For greater angles, the torque 
per current can be several factors higher, a result that 
seems counterintuitive. In fact, as the barrier is made 



3 




15 20 5 10 15 20 

layer index in free Fe layer 

FIG. 3: (Color online) Layer resolved torque per current 
acting on the free layer at S = 90 °. Left: Free electron 
model using a simple rectangular barrier (B) with an imag- 
inary wavevector oi kg = kg = 4 i nm~^ for majority (t) 
and minority spins (J.) within the barrier. Left (L) and right 
(R) of the barrier we use the half-metallic case present in 
Fe at the F point. That is, the majority wavevector is real 
nm , whereas the minority wavevector is 
imaginary k\^ = k^ = 5 i nm~^. The thickness of the barrier 
is 80 nm. Right: For comparison the ab initio result is shown 
integrated over the whole Brillouin zone, for a free layer with 
20 monolayers. 



thicker, there is no intrinsic limit to how large the quan- 
tity can get. This behavior results from the contribution 
of electrons that penetrate the barrier but reflect. They 
contribute to the torkance but not to the conductivity. 
Consider a free layer with a magnetization at a relative 
orientation of either 45° or 135° with respect to that 
of the fixed layer. In both cases, the majority electrons 
that tunnel from the fixed layer have the same transverse 
component and will exert roughly the same torque. The 
electrons tunneling into the free layer with a magnetiza- 
tion at 45° are approximately 85 % majority and 15 % 
minority in the free layer while the electrons tunneling 
into the free layer with a magnetization at 135° are ap- 
proximately 15 % majority and 85 % minority in the free 
layer. The tunneling current will be much greater in the 
former case than in the latter so that the ratio of the 
torque to the tunneling current will be greater when the 
layers are at 135°. 

The weak infiuence of free layer thickness on dr^^/dl 
in Fe/MgO/Fe and the saturation for more than three 
monolayers of Fe implies that the torque is restricted to 
the layers next to the interface even though there is only 
minimal dephasing. To understand this behavior, the 
right panel of Fig. [3] shows our ab initio results for the 
layer resolved torque. We analyze this result by compar- 
ison to a simple free electron model in the left panel of 
Fig.H 

In the simple model, we use the wave vectors of the Ai 
band in Fe at the Fermi level to model the half metallic 
nature at the F point. In a half metal, the components 
of the non-equilibrium magnetization perpendicular to 
the magnetization and hence the torque arise from co- 
herent interference between the majority and minority 
spin components. Since the minority state is evanescent. 



the torque must decay. After a short distance only the 
component of the spin current pointing along Mn is left. 
The decay of the torque, which results from the half- 
metallic character of the ferromagnetic electrodes, is even 
faster than in spin- valve systems where the decay follows 
a power law resulting from dephasing [3]. The agree- 
ment between the ab initio results and the simple model 
emphasizes that the states around the Brillouin zone cen- 
ter dominate the torque within the free layer. 

During this decay, the electrons precess around Mn 
within the zy plane leading to decaying oscillations of the 
components of the magnetic moments of the conduction 
electrons along the z and y directions. These oscillations 
in the magnetic moments are phase-shifted by 90 ° with 
respect to each other [2^ . They lead to the oscillations in 
the torque components having the corresponding phase 
shift as discussed in Fig. [2l Within a given layer, the 
sizes of r|| and t± are comparable on average. The total 
torques are determined by the starting phase and the 
strength of the decay. In addition. Fig. [3] only shows 
the torque of right-going electrons. The contributions of 
left-going holes have to be added in order to obtain the 
torque in Eq. [TJ 

There are additional contributions to the torque that 
are important for thinner barriers but are negligible for 6 
monolayer and thicker barriers. In Fig. [21 there are oscil- 
lations that are barely visible, which originate in states 
away from the Brillouin zone center for which Fe is not 
half-metallic. These oscillations are long ranged due to 
reduced dephasing and could be important for barriers 
that are thinner than those typically fabricated. 

Our calculations of the bias dependence of the out- 
of-plane component of the torkance address a topic of 
recent experimental interest. Spin-transfer-driven ferro- 
magnetic resonance (ST-FMR) quantitatively measures 
the magnitude and direction of the spin transfer torque 
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in tunnel junctions [271 . |28| . I29j . Tulapurkar et al. 
measured a linear dependence of t±_ on the applied bias 
for small voltages (non-zero dT±/dV) whereas Sankey et 
al. [s^l measured t± to be linearly independent of V. 
This latter result is consistent with theoretical argu- 
that hold in special cases 
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32| that hold in special cases. For junctions 
that are symmetric about the center of the barrier and in 
which electrons leaving the ferromagnetic layers into the 
leads are aligned with the magnetizations, the left going 
electrons cause the same out-of-plane torque in the free 
layer as the right going ones, that is Tj^{Ep) = t^{Ef)- 
When the outgoing electrons are aligned with the magne- 
tizations, the total torque on both layers in a symmetric 
junction must be in-plane. If this condition holds, Eq. [1] 
gives dT±/dV = and there is no linear dependence on 
the applied voltage for small biases. However, for asym- 
metric junctions these arguments do not apply. 

Figure [2] shows that it is possible to get a significant 
out-of-plane component of the torkance in tunnel junc- 
tions that are close to ideal. However, thickness fluctua- 
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FIG. 4: (Color online) Thickness fluctuations. Left: Sketch of 
assumed macroscopic fluctuation of Fe layer thickness. Fixed 
layer consists of equal parts 19 monolayers and 20 monolayers 
of Fe. The free layer consists of 50 % n, 25 % (n + 1), and 
25 % (n— 1) Fe monolayers. Right: Torkance as a function of 
the averaged free layer thickness n assuming that each of the 
six possible junctions can be described separately. The total 
torkance is calculated by conducting the junctions in parallel. 
Note that the scale is identical to the scale in the bottom left 
panel of Fig. [2] 



tions in a real junctions lead to an averaged torkance. To 
simulate such fluctuations we assume the structure shown 
in the left panel of Fig. H) The right panel shows the 
corresponding averaged torkance. Comparing to Fig. [2] 
shows two consequences of the thickness fluctuations. 
First, the oscillation of the in-plane component vanishes 
and second, the out-of-plane component itself vanishes. 
This behavior is consistent with the experimental re- 
sults in Ref. [s^. Therefore, to actually measure a non- 
vanishing out-of-plane component of the torkance or an 
oscillating behavior of the in-plane component requires 
samples that are close to ideal. For most samples the 
out-of-plane component will vanish. 

In conclusion, we show that the spin-transfer torque 
in Fe/MgO/Fe tunnel junctions behaves very similarly 
to the behavior found in all-metallic devices. The spin 
transfer torque is largely localized to the interfaces and 
largely in the plane defined by the two magnetizations. 
The dominant contribution to the tunneling and the 
torque comes from states around the Brillouin zone cen- 
ter where Fe is a half-metal with respect to the Ai states. 
This half-metallic behavior leads to an exponential decay 
of the torque within the ferromagnetic layer even without 
dephasing. For a perfect sample we expect some small 
out-of-plane component and an oscillation of the in-plane 
component of the torkance as a function of the free layer 
thickness. However, small fluctuations of the thickness 
will average out this component. 
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